Effects of Disorder and Interactions in the Quantum Hall Ferromagnet 
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This work treats the effects of disorder and interactions in a quantum Hall ferromagnet, which 
is realized in a two-dimensional electron gas (2DEG) in a perpendicular magnetic field at Landau 
level filling factor v = 1. We study the problem by projecting the original fermionic Hamiltonian 
into magnon states, which behave as bosons in the vicinity of the ferromagnetic ground state. The 
approach permits the reformulation of a strongly interacting model into a non-interacting one. The 
latter is a non-perturbative scheme that consists in treating the two-particle neutral excitations of 
the electron system as a bosonic single-particle. Indeed, the employment of bosonization facilitates 
the inclusion of disorder in the study of the system. It has been shown previously that disorder 
may drive a quantum phase transition in the Hall ferromagnet. However, such studies have been 
either carried out in the framework of nonlinear sigma model, as an effective low-energy theory, 
or included the long-range Coulomb interaction in a quantum description only up to the Hartree- 
Fock level. Here, we establish the occurrence of a disorder-driven quantum phase transition from a 
ferromagnetic 2DEG to a spin glass phase by taking into account interactions between electrons up 
to the random phase approximation level in a fully quantum description. 

PACS numbers: 



I. INTRODUCTION 

The simultaneous treatment of disorder and interac- 
tions in strongly correlated electron systems has always 
formed a knotty challenge; this is because of the dearth of 
manageable analytical techniques that can deal with dis- 
order and interactions at the same time.— ^ The strongly 
correlated system of interest in this work is the two- 
dimensional electron gas (2DEG) in a perpendicular mag- 
netic field at Landau level filling factor v = 1, whose 
ground state is commonly known as the quantum Hall 
ferromagnet. 

The quantum Hall ferromagnet is the spin-polarized 
ground state of the 2DEG at v = 1 in which all electrons 
completely fill the lowest Landau level with spin up po- 
larization. Such configuration minimizes the Coulomb 
energy for fermionic systems. In general, it is a competi- 
tion between kinetic and Coulomb energies, which deter- 
mines the ground state. In the case of the quantum Hall 
ferromagnet having v = 1 the kinetic energy is frozen and 
does not change with spin flip, thus, the ground state is 
ferromagnetic, even with zero Zeeman splitting. Typi- 
cally, the Zeeman splitting in the GaAs heterojunctions 
turns out to be roughly 70 times smaller than the spacing 
between Landau levels and an order of magnitude smaller 
than the Coulomb energy per particle. 

The neutral elementary excitations are spin wave ex- 
citations, also called magnons. The spin waves can be 
described by the action of the spin lowering operator 
Sq, projected to the lowest Landau level, on the ferro- 
magnetic ground state. It turns out, that the projected 
operator creates an exact excited eigenstate of the Hamil- 
tonian. In the regime of low momenta, the magnon's dis- 
persion is quadratic and the coefficient of the quadratic 
term represents a phenomenological constant known as 



the spin stiffness. The spin stiffness provides a measure 
of the free-energy increment associated with twisting the 
direction of the spins. A significant spin stiffness indi- 
cates that the system lies in the ferromagnetic phase, 
while a paramagnetic state corresponds to a vanishing 
spin stiffness. The spin wave dispersion at very large 
momenta saturates at a constant value given by the sum 
of the Coulomb and Zeeman energies. Thus, at large mo- 
menta, the value corresponds to the energy of separate 
quasiparticle and quasihole excitations. 

One approach that has successfully dealt with strongly 
correlated electron systems is the so-called bosonization 
procedure. Bosonization is a non-perturbative approx- 
imation scheme that essentially treats the electron-hole 
excitation, known as exciton, as a bosonic single-particle; 
consequently, a fermionic Hamiltonian can be recast into 
a bosonic one. In 1950, Tomonaga revealed, in a ground- 
breaking paper,— that the application of the bosoniza- 
tion formalism to a one-dimensional electron gas (1DEG) 
yielded an exactly-solvable Hamiltonian. The reason is 
that the electron and the hole propagate with nearly the 
same group velocity in the low-energy region. However, 
that is not the case in two dimensions. At a given momen- 
tum k, the particle-hole pair excitation holds a continu- 
ous range of energies. Therefore, it is less straightforward 
to construct a coherently propagating bosonic entity in 
two dimensions. 

The first attempt to extend the bosonization procedure 
for higher dimensions was done by Luther— and then re- 
vised by Haldane£ Castro Neto and Fradkinji as well 
as Houghton and Marston,— & developed a bosonization 
technique for a Fermi liquid in any number of dimen- 
sions. As regards the interacting 2DEG subject to an 
external perpendicular magnetic field, Westfahl Jr. et 
a/.— constructed a formalism that treated the elemen- 
tary neutral excitations of the system, the magnons, in a 
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bosonic framework such that the fermionic Hamiltonian 
of the system was transmuted into a quadratic bosonic 
Hamiltonian. The drawback is that this method is valid 
in the limit of weak magnetic fields, which amounts to 
large Landau level filling factors v. 

Doretto et alM extended the methodology of West- 
fahl Jr. et alr£ to the case of the 2DEG at v = 1 
(i.e. for a very strong magnetic field). Given that the 
system is now restricted to one Landau level, the task 
greatly simplifies, since the Landau level quantum de- 
gree of freedom can then be disregarded. Projecting the 
original fermionic interacting Hamiltonian of the system 
into the lowest Landau level, which is completely filled 
(v = 1), allows one to expand it in magnon states.— It 
then turns out remarkably that the dispersion relation 
of the free magnons coincides with the result derived 
by Kallin and HalperinAi within the fermionic descrip- 
tion at the random phase approximation (RPA) level and 
the quartic interacting part of the magnon Hamiltonian 
might be related to the skyrmion-antiskyrmion neutral 
excitations of the Hall ferromagnet.— Moreover, in the 
vicinity of the ground state, without magnon-magnon in- 
teractions, magnons behave like bosons. This allows to 
treat magnons approximately as bosons in the so-called 
single-mode approximation. ~ 

Here, we intend to calculate a quantum phase tran- 
sition in the quantum Hall ferromagnet driven by dis- 
order, accounting for the Coulomb interactions between 
electrons. We will use the bosonization technique allied 
to the usual self-consistent Born-approximation for the 
disorder averaging procedure. 

Before presenting the results obtained in this paper, 
it is worth getting acquainted with the current status 
of research related to the field. To begin with, Green 1 
propounded that the vanishing of the renormalized spin 
stiffness at a threshold value of the disorder strength sig- 
nifies the occurrence of a depolarization transition from 
the ferromagnetic phase to a paramagnetic one. His find- 
ing is based upon a previous result established by Fogler 
and Shklovskii, 15 who proffered the same idea in the case 
of higher Landau levels. Green established this proposi- 
tion in the framework of non-linear sigma model, used 
as an effective low-energy theory in the regime of weak 
disorder. The other quantity that Green computed is the 
disorder contribution to the optical conductivity, which 
he found to be unmeasurably small. Finally, Green es- 
tablished that the quantization of the Hall conductivity 
is not affected by the presence of weak disorder in the 
system.— 

Another work was carried out by Sinova, MacDonald 
and GirvinyS. who established the occurrence of a phase 
transition from the paramagnetic state to the partially- 
polarized ferromagnetic one and then finally to the fully- 
polarized ferromagnetic one as the interaction strength 
increases relative to the disorder strength. They deter- 
mined this result by computing the average value of the 
spin polarization as a function of the interaction strength 
relative to the disorder strength. Sinova et al£ did con- 



sider Coulomb interactions within the framework of the 
Hartree-Fock approximation. Moreover, the transition 
from the paramagnetic phase to the ferromagnetic one 
was found to take place when the Coulomb energy scale is 
about twice as large as the Landau-level-broadening dis- 
order energy scale. As a final point, the authors inferred 
that no phase transition can take place in the strong dis- 
order limit. 

The last germane paper was published by Rapsch, Lee 
and Chalker.— They established the occurrence of a phase 
transition from the ferromagnetic state to the so-called 
spin glass phase. This result was obtained by calculat- 
ing the magnetization, the magnetic susceptibility and 
the spin stiffness as functions of the disorder strength. 
They assumed the disorder potential to be Gaussian dis- 
tributed and described the system in terms of a semiclas- 
sical spin model. In their model, they took into account 
Coulomb interactions within the Hartree-Fock approxi- 
mation but modelled them as being short-ranged. Like 
Greenyi Rapsch et aly computed the disorder contribu- 
tion to the optical conductivity and found as well that 
it is undetectable. Finally, they calculated the dielectric 
susceptibility of both the partially-polarized ferromag- 
netic phase and the spin glass one and they concluded 
that both regimes display an insulating behavior at low 
momenta and a metallic behavior at large momenta. 

Let us now put our work in perspective. Our objec- 
tive is to establish the behavior of the renormalized spin 
stiffness as a function of the disorder strength in order 
to ascertain a potential quantum phase transition driven 
by disorder to a non-ferromagnetic state. Indeed, if the 
spin stiffness vanishes for a critical value of the disorder 
strength, then this signals an instability in the ferromag- 
netic phase.— On the other hand, the appearance of an 
imaginary component of the spin stiffness, which might 
be interpreted as a spin wave damping , 16 ' 17 at a certain 
disorder strength, might indicate the appearance of lo- 
calized spin waves and a spin-glass phase transition. An- 
other important characteristic is the Pauli susceptibility, 
which diverges at the point of the phase transition from 
a non-ferromagnetic to a ferromagnetic state, indicating 
spontaneous magnetization. We consider a fully quantum 
model, include a short-range weak disorder potential up 
to the 2 nd order Born approximation and treat the true 
long-range Coulomb interactions up to the RPA level. 

The method that we employ consists of five steps. First 
of all, a bosonized expression of the total Hamiltonian, 
which includes a contribution from disorder, is sought 
for. The dispersion relation of the free bosons corre- 
sponds to the one computed by MacDonald et alr& and 
more explicitly by Doretto et al.^ which entails interac- 
tions between electrons up to the RPA level. The second 
step consists in obtaining the full Green's function, and 
precisely its disorder self-average. In our case, because 
the impurities are randomly distributed throughout the 
system, the disorder self-average can also be taken by 
averaging over the impurity positions. The third stage 
is then to determine the self-energy of that disorder self- 
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averaged Green's function through the use of the Dyson's 
equation. The self-energy is determined in the low- 
impurity density and weak disorder scattering approxi- 
mations. As a result, the self-energy corresponds to a 
single diagram with one propagator line and two disorder 
potential lines. The propagator line is evaluated within 
two further possible approximations: the bare approxi- 
mation, which consists in using the bare bosonic propa- 
gator, and the self-consistent approximation, which uses 
instead the full disorder self-averaged Green's function. 
One must bear in mind that both propagators take into 
account interactions between electrons up to the RPA 
level. Furthermore, the bare approximation is first taken 
in the long wavelength limit, which keeps the lowest order 
terms in momenta, and then in the general case, where all 
the momenta terms are taken into account. The fourth 
step consists in obtaining the renormalized dispersion in 
these approximations: bare and self-consistent approxi- 
mations. The final stage is then to determine the spin 
stiffness in the approximations by taking the coefficient 
of the quadratic term in the renormalized dispersion. It 
is found that a naive extrapolation of the bare approxi- 
mation to the regime of finite disorder strength predicts 
vanishing of the renormalized spin stiffness at a certain 
disorder strength u p , indicating a paramagnetic phase 
transition. A more realistic self-consistent approxima- 
tion, however, predicts even faster decrease of the renor- 
malized spin stiffness with growing disorder strength up 
to a certain critical value u c of the disorder. At this 
point, the renormalized spin stiffness drastically changes 
its behavior: it becomes nonanalytic, acquires an imagi- 
nary part, and the real part saturates at a certain positive 
value without reaching zero. Such nonanalytic behavior 
cannot be accessed by any finite number of perturbative 
corrections. In addition, our calculations show a strong 
indication that the Pauli susceptibility also diverges at 
the same critical point u c , suggesting a phase transition, 
presumably to a spin glass phase. 

The outline of this paper is the following: in Section 
II we present the model and in Section III we derive the 
expression for the self-energy. Then, we first solve the 
problem using the bare Green's function in Section IV. 
We present our numerical and analytical results for the 
self-consistent solution of the Dyson equations in Section 
V and draw our conclusions in Section VI. 
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Wimp - [dr J2 u (* ~ Xi)* f (r)*(r). 

i—l 

Here, \&t(r) and \l/(r) are, respectively, the fermionic cre- 
ation and annihilation operators in coordinate space, m* 
denotes the effective mass of the electron, A is the vector 
potential, g* stands for the effective Lande g-factor and 
fiB is the Bohr magneton. In addition, V(|r|) = e 2 /(e|r|) 
denotes the Coulomb potential, with e being the dielec- 
tric constant of the host semiconductor, and U stands 
for the impurity potential, with Xi being the random 
position of an impurity. 

The first step consists in obtaining a 2 nd quantized 
version of the magnon Hamiltonian of the system. In our 
model we consider only single magnon processes, which 
allow us to use a bosonic description. It was shown in 
Ref. that the bosonized Hamiltonian of the system 
in the absence of disorder is (neglecting a constant term) 



(1) 



where 6 q and 6 q are, respectively, the bosonic creation 
and annihilation operators in q space and the bosonic 
dispersion relation is given by 



Wq = g + (-B 



1 _ e -l^l 2 /4 Jn 



N 2 



(2) 



Here, cb = y / n/2(e 2 /el) stands for the Coulomb energy 
scale it being the magnetic length), Iq denotes the mod- 
ified Bessel function of the first kind, and g = g* [IbB. It 
must be stressed that although the interaction between 
magnons is omitted from the discussion, the Coulomb in- 
teraction between electrons up to RPA level is taken into 
account by the bosonic dispersion relation <x>q i n i 13 

We now focus on the impurity part of the Hamilto- 
nian. We begin with the fermionic expression of the 2 nd 
quantized impurity Hamiltonian, 



II. THE MODEL 



Wimp = J2 U fl pH 



(3) 



The 2DEG in the presence of both a perpendicular 
magnetic field (B = Bz) at v = 1 and disorder is de- 
scribed by the fermionic Hamiltonian 7i = Ho + 7Yi mp , 



Here, <? q denotes the Fourier transformed density func- 
tion X^T P ^( x — Xj) for the impurities and a q and a q 
are, respectively, the fermionic creation and annihilation 
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operators in q space. In order to obtain the bosonic form 
of the above, the Fourier-transformed electronic density 
operator must be used. It is given by 



P(q) = / dr e _iq r * t (r)«'(r). 



(4) 



The electronic field operators are related to the single- 
electron operators by 

g— ipr e ipr 

*( r ) = z3^=r a p and ^( r ) = l2-77 a l 

where ^4 is the area of the system. Substituting the above 
back into Eq. (|4| gives 



p(q) 



dr e- Iq - r ^ 



D i(p'-p)-r 

' 



t _ \ " t 

p' a P ~~ / j a p+q Q P- 



(5) 



Then, substituting Eq. back into Eq. (J3j) yields 

Wimp = 53 ^ qp(q)- ( 6 ) 
q 

The bosonized version of the electron density operator 
reads^i 

P(q) - S^ N, + 2l e-l^l 2 / 4 53sin (5^2) & q+p&p , (7) 

p 

where jV^ = A/(2irl 2 ) is the Landau level degeneracy 
and qAp = £ 2 z ■ (q x p). The disorder Hamiltonian then 
becomes 



Wimp = E q ^(q)3q [V^<A + 2ie^q| 2 /4 

xE p sin(^)& q+p & P 



(8) 



The constant term <5 qi o^</> is now omitted since the quan- 
tity of interest is the Green's function. 

The bosonized impurity Hamiltonian is then finally 
written as 



Wimp = 53 C/ (q)5q/(q i P) & q+p&P> 



where 



/(q,p) = 2 ie -l^/ 4sin (SAP) 



Labelling 



t/(q)/(q,p) = £/ e (q,p), 



(9) 



(10) 



(11) 



the full bosonized Hamiltonian of the quantum Hall fer- 
romagnet in the presence of impurities is then expressed 
as 

n = 53 ^ q 6 q 6 q + 53 C/ e (q, P ) 5q & q+P V (12) 



Let us now say a few words on the dimensions of the 
disorder potential. There are two sources of disorder 
present in the system: impurities positioned at a cer- 
tain distance away from the 2DEG and impurities present 
in the 2DEG. In the case of GaAs heterostructuresji^ 
most of the disorder potential is spawned by the Coulomb 
interaction between the electrons and the impurities lo- 
cated away from the 2DEG. These impurities correspond 
to ionized donor atoms situated in the n-type region, 
which itself is detached from the 2DEG by an insulat- 
ing layer of thickness d ~ 1000A ^> I. In the present 
calculations, the disorder potential will be taken as an 
effective two-dimensional potential. 

Having obtained the bosonized Hamiltonian in the 
presence of impurities, one is now able to determine the 
expression for the self-energy. 



III. DERIVATION OF THE SELF-ENERGY 

In the same spirit as Ref [20], one first looks for the 
Green's function, 

G(p', P ;i) = -i(0|T[6 p (i)6 p/ (0)]|0). (13) 

Here, |0) stands for the bosonic vacuum state, which is 
none other than the quantum Hall ferromagnet: i.e. |0) = 
\QHF) = Ult^cl^F- Thus, one has 

z^G(p', P ;t)= ( 5(t) ( 5 p , p ,-*(0|T 1 [[6 P (t),W]6 p ,(0)] |0), 

(14) 

where T is the time ordering operator. Now, using 
Eq. (fl2"ll . one easily finds that 



[b p (t),H] =^ p 6 p (t)+53[/ e (q !P -q). gq 6 p _ q (i), (15) 
q 

such that one obtains for the second term in Eq. (fl4| 
-i(0\T[[b p (t),H]bl,(0)] |0)=WpG(p',p;t) 

+ E q £ /e (q:P-q).9qG(p',P-q;i)- (ie) 

Hence, the equation of motion of G(p', p; t) is written as 

d 



ifa-up) G(p', p; t) = S p , v >S(t) 



53 c/e ( q ' p ~ q)5q G (p'> p - q; *)• 



(17) 



The zero-order approximation to the solution of 
Eq. (TU yields 



G°(p',P;*) = W<?°(p,i), 



(18) 



where G°(p,t) stands for the bare bosonic Green's func- 
tion. We now look for the expression for G°(p,i). 

Firstly, one needs to find the Heisenberg bosonic op- 
erator in the absence of the disorder potential. Starting 
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<G(p)>. 



p 



p P 




p P 




U (q-p,p) \ / y e (p . q , q , 
N. 



- u (q'-h.q) 

V | ✓ 

u" (q-p,p) \ ! »' 1/ (p-q',q') 

V 

N. 



Figure 1: Diagrammatic expansion of the disorder averaged Green's function. 



with id t b p (t) — [b p (t),Ho] — ui p b p (t), one then obtains 
b p (t) = & p e~ lw p*. Therefore, for the case t > 0, the free 
Green's function is 

G°(p,t) = -<(0|& t ,(t)6t|o) = _i e -*«P*<o|6 p 6t|o) 



Here, the bare Green's function reads 



G (p,w) 



1 



uj — uj p + it] 



(21) 



whereas for t < 0, it turns out to be 

G°(p,t) = -*<0|6pC*)6t|0) = -ie-^*(0|6t fep |0) = 0. 

This solution is indeed identical to the electronic one. 

Now, the cynosure is on the generic solution of the 
differential equation (fTF)) . By coupling the latter with 
the boundary equation (fT8|l yields the integral equation 



/oc 
df G°(p,t-t') 
-OO 



E C/6 ( c l' P ~ q)5qG(p', P - q; i). 



By Fourier-transforming the time in Eq. (fT9|) to fre- 
quency and shifting q — » p — q one finds 

G(p', p; w) = <5 P , P 'G°(p, lo) + G°(p, w) 

X^C/ e (p-q,q)g P - q G(p',q;^). ( 2 °) 



where 77 — > + and u> p is given by Eq. $2$. In the same 
way as for fermions, the solution of Eq. (|20|) is obtained 
by iteration. One gets the so-called Born series: 



g(p'p) = E gW (p'p). 



(22) 



(19) where G°(p',p) = 5 p ^G°(p) and for n> 1, 



G {n) (p', p) = G°(p) U e (p ~ q, ^9 V ^G^ (p', q). 



Expansion of Eq. ([22)1 then yields 



G(p', p) = ^p,p'G°(p') + G°(p')C/ e (p - p', p').9p-p'G°(p) + ]T G°(p')C/ e (q - p', p').9 q -p'G (q)(7 e (p - q, q).9 p - q G (p) 

q 

+ J2 G°(p')U e (q - p', p').9q- P 'G°(q)L/ e (q' - q, q). 9q '- q G (q')£/ e (p - q', q')Sp- q 'G (p) + . . . , 
q,q' 

Due to disorder self-averaging in the limit of very large number of impurities N- lmp — > 00, with constant density 
"imp = const., the full bosonic one-particle Green's function approaches its average value 

((G(p', p) - <G(p', p)) im p) 2 ) im p - 0, (23) 

which is 

(G( P ',p)) imp = WG°( P ') + <.g P -p')impG°(p')f/ e (p - p',p')G°(p) 

+ E(-9q-P'.9p-q>im P G°(p')t/ e (q - p', p')G°(q)C/ e (p - q, q)G°(p) 
q 

+ E^q-P'3q'-q5p-q')impG (p')C/ e (q- p',p')G (q)(7 e (q' - q,q)G°(q');y e (p - q',q')G°(p) + . . . 
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In the thermodynamic limit A — > oo: 

(^q)imp -^imp^q,0i 
(ffqSp)imp = ^hnp^P.O^q.O + Mmp#q+p,0, 
(3q'Sq#p)imp = -^knp^q'.O^q.O^p.O + -^hnp (<Wq,0<5q' ,0 + <5q+q',0^p,0 + ^p+q',0<5 q ,o) + -Wimp<5q'+q+p,0 • (24) 



Moreover, one has 



C/ e (0,p) = f7(0)/(0,p) = C/(0)2ie-l £ (°)l 2 / 4 sin ^ 



OAp 



0. 



(25) 



Substituting Eqs. ([24]) and Eq. lj25|) into the expression for (G(p))i mp shows that the translational invariance is recov- 
ered after the averaging (G(p',p))i mp = (G(p))i mp 5 p / )P , where 



(G(p)) imp = G°(p) + N imp G°(p)C/ e (q - p, p)G°(q)f/ e (p - q, q)G°(p) 

q 

+ N imp ]T G°( P )^(q " P, p)G°(q)C/ e (q' - q, q)G°(q')^(P - q', q')G°(p) 



(26) 



Therefore, there is no 1 st order Born scattering contribution to the bosonic self-energy. Moreover, it is possible to 
show that all odd order contributions to the self-energy vani sh (see Appendix A). 



Bp) 



q 



U e (q-p,p) \ / U e (p-q,q) 
N 

imp 

Figure 2: Self-energy in the low-density weak scattering ap- 
proximation. 



low-density weak scattering approximation will be used 
throughout the calculations. Low density means that the 
number of disorder atoms present in the system is taken 
to be much lower than the number of electrons, while the 
weak-scattering approximation signifies that the scatter- 
ing potential induced by a given impurity atom is weak, 
such that only the first and second-order Born scatterings 
are accounted for. The problem then reduces to solving 
the diagrammatic expression shown in Fig. [51 



This result is expressed diagrammatically in Fig. [1] It 
was shown2& that the disorder averaged Green's function 
can also be expressed as 



(G(p)> 



1 



imp 



(27) 



*-Wp- S(p,w) 
Hence, the self-energy must now be computed. The 
I 



The self-energy can be evaluated in two different man- 
ners: 1) the bare approximation that uses the bare 
propagator G° and 2) the self-consistent approximation 
that uses the full disorder self-averaged Green's function 
(G)i m p. Therefore, for generality we will use the propa- 
gator G(q, oj), which is going to be specified further for 
each particular case. This yields algebraically 



E(p, u) = N imp U e (q - p, p)G(q, uj)U c (p - q, q) 
q 

= N imp ^2 u (l ~ P)/(q - P> P)G(q, u)U(p- q)/(p q, q), 



(28) 



where Eq. (fTTj) was substituted in the second line. 

In this work, the impurity potential is assumed to be short-range, i.e. U(q) — constant. An uniform potential 
in momentum space is attained from a delta function interaction in real space, U(r) — U6(r) (such that E/(q) = 
(l/A) J dre lqr [/<5(r) = U/A). Thus, this model assumes that the bosons (and therefore the electrons) collide directly 
with the impurity 'atoms' that constitute the effective disorder potential; in reality, most of the impurities are located 
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away from the 2DEG. Thus, one has 



E(p,w) = N imp ^2 



lie 



efq-p)l 2 /4 dn f (q-p)Ap 



G(q,w) 



2ie -K(p-q)| 2 /4 



sin 



(p - q) A q 



4iVi 



imp 



.4 



-|*(q-p)| 2 /2. 



sin (__jG(q,w). 



(29) 



One then expands the argument of the exponential: 



e -K(q-p)| 2 /2 = e -Kq| 2 /2 e -|^p| 2 /2 e f 2 q-p = g" {l<lf /2g- (^p) 2 /2^ 2 gp cos c 



(30) 



Here, denotes the angle between vectors q and p. Furthermore, the summation is transmuted into an integration 
through the use of the formula, 



A '" x 
4^2 



dq q / dO. 

II JO 



(31) 



The angle 9 is taken arbitrarily on the plane containing the vector q, therefore, one is free to set = <ft. The sine 
squared term in Eq. ([29|) can be re-written as sin 2 (q A p/2) = [1 — cos(q A p)]/2 = [1 — cos(^ 2 gpsin</>)]/2. We also 
assume rotation invariance of the Green's function G(q, lj) = G{q,us). Hence, the self-energy is also rotation invariant 
and can be expressed as 



E(p,cj) = 4n imp [/ 2 



^1 ge -(^) 2 /2 p -(^) 2 /2 



2tt 



G(q,oj) e i2 ™ cos *l[l-cos(£ 2 qp sin^)], 



(32) 



where n- unp = N- imp /A stands for the impurity density. After a straightforward calculation (see Appendix B), we find 

S(p, W ) = 4n imp [/ 2 ^°° g g e -W 3 /2 e -W 2 /2 G(9)W )l [7o( ^^ ) _ !]. (33) 

Rescaling the momenta by q, p — > q/£, p/£ simplifies the self-energy to 

poo 

nP,u) = - A e 2 B e-v 2 ' 2 J dqqe-i 2 / 2 lI (qp)-l]G(q 7 u;) 7 (34) 

where the various pre-factors, including the disorder potential strength and the impurity density, can be re-grouped 
into a single convenient parameter: 



4t7- U 2 



(35) 



which will be dubbed the disorder strength. Thus, u is a dimensionless parameter that measures the disorder interaction 
strength relative to the Coulomb interaction, u ~ (Edi s / E cou \) 2 . The above self-energy expression will be evaluated 
in two different ways: i) first order corrections in u and ii) self- consistently. 



A. Bare Approximation 

In the bare approximation the self-energy (|34| becomes 



" 2 ^l 2 I dqqe-f/ 2 [Io(qp)-l]G°(q,Lj). 



After substituting Eq. 11211) into the above, we obtain 



dq qe 



V/ 2 hjqv) - 1 

U) — UJq + IT) 



(36) 



(37) 
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Figure 3: (Color online) (a) Renormalized dispersion in the bare full k approximation (blue), in units of the Coulomb energy 
e 2 /(e£), contrasted with the one in the bare long wavelength approximation (red), both as functions of the momentum \p£\ and at 
u — 0.1. (b) Renormalized spin stiffness in the bare full k approximation (blue) and in the bare long wavelength approximation 
(red). Notice that using the bare Green's function Go we find a transition from a ferromagnetic to a paramagnetic phase, 
whereas using G in the self-consistent approximation we find a transition into a spin glass phase (see next section). 



Making use of the identity (for r\ — > H 



= V in8(x), 

x + vr\ x 



we find the real and imaginary parts of the self-energy 



ReS(p,w) 



Im£(p, ui) 



"P 2 /2 



-9 2 /2 



h(qp) 



VI dq qe 

LJ-LJ q 



V/2 



dq qe 



V/2 



[I (qp) - 1]tc6(u - U q ). 



(38) 

(39) 
(40) 



The above equations can be evaluated analytically in the long wavelength approximation, which is done in Appendix 
C. Here, one uses the complete bosonic dispersion relation given by Eq. |(2]). As a result, one can only solve the 
imaginary self-energy numerically; that task is not performed here. We concentrate, instead, on the real part. 

The renormalized energy of the bosons (including the disorder contribution) is obtained by looking at the poles 
of the full disorder self-averaged Green's function in Eq. f27j) . u> — u> p — ReS(p, uS) = 0, such that the renormalized 
dispersion relation is determined from Eq. f39j) : 



uj = g + e B 



1 



- p2/2 V / dq qe 
Jo 



9 2 /2. 



Io(qp) - 1 



9 + e B l-e-« 2 /4/ * 



(41) 



The corresponding plot is illustrated on Fig.[3](a). One can notice that at not too large momenta (i.e. near \p£\ = 1) 
there exists already a substantial difference between the bare (long wavelength) and bare (full k) approximations. 

Now, the renormalized spin stiffness is sought for. For the sake of convenience, one begins by introducing the 
variables ui,g = ui / cb, g / 'eg and re-writing Eq. f4Tj) as 



g = i-e-^%(^) + y^vl°°dq 



qe 



-q 2 /2_ 



h{qp) - 1 



(u) - g) ~ 



1 - e-? 2 / 4 I * 



(42) 



One then expands the above in powers of p, 

,2 



1 



El 

2 



64 



V I dq qe' 

o 



V/2. 



(2 2 



- 9) - 



(43) 
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and one takes only the p 2 terms, 



P 2 


u p 2 


4 


~ IT 


P 2 


[i-- 

4 


4 



g 3 e -? 2 /2 



1 - e-9 2 /4/ ( 



■•(*) 



(44) 



(5.72) 



The renormalized spin stiffness then reads 



Ps 



T' 



1 - 1.43m). 



(45) 



Eq. (|45|) is the main result of this section. The above expression was derived in the bare approximation, which takes 
into account only the lowest order corrections in u. Such assumption is only true in the realm of weak-disorder 
scattering.— It can be seen that the renormalized spin stiffness decreases linearly in this approximation. A naive 
extrapolation of this dependence to the region of finite and strong disorder strength shows that there is a certain 
value Up — 0.7, for which the renormalized spin stiffness vanishes (in the long wavelength approximation u p — 1), see 
Fig. [3jb) . Green- explains that a vanishing renormalized spin stiffness at a threshold disorder strength means that the 
2DEG at v = 1 undergoes a quantum phase transition from a ferromagnetic state to a paramagnetic one. Thus, one 
can infer that the quantum Hall ferromagnet undergoes a disorder-driven quantum phase transition to a paramagnetic 
state at critical disorder strength u p = 0.7. It is also interesting to remark that Green established this general finding 
in the domain of the weak disorder limit (though in the context of a different model) . The results obtained in this 
section cannot be directly compared quantitatively with those of Green^, Sinova et al£ and Rapsch et al? In addition 
to the fact that the model used in the studies of Green is different, he does not complement his proposition on the 
vanishing of the renormalized spin stiffness with some quantitative results. Sinova et al. 2 use a disparate variable in 
the ratio of the interaction strength to the Landau-level broadening disorder energy scale. Finally, Rapsch et a/.— 
perform their numerical calculations on a semiclassical spin model. 

In the next section we evaluate the self-energy using the so-called self-consistent approximation and show that the 
renormalized spin stiffness drastically changes its behavior, which leads to completely different conclusions about the 
phase transition. 



IV. SELF-CONSISTENT APPROXIMATION 



The self-consistent approximation means that the self-energy is evaluated with the total disorder averaged Green's 
function lj2"7j) instead of the bare one. Therefore, one has (see Eq. ((34J) ) 



POO 

Y> u {pM = -/B^ 12 J dqqe-« 2 / 2 [I (qp)-l}(G u (q,oj)) inip , (46) 
Now, by referring to the computations carried out in the previous section and substituting Eq. (|27ll . one gets 



Using that 



one has 



V ( \ U 2 -p 2 /2 f i -q 2 /2 A)(<?P)-1 

E u (p,w) = -e B e p ' I dqqe q/ — -. 47 

4 J oj - cj 9 - £„(<?, w) 



m?p)-i=eS. ( 48 ) 

n— 1 ^ ' 



11 2, ^ T) 2n r°° „2n+l - 9 2 /2 

E u (p, u) = ^e 2 B e-P / 2 Y / dq — ? ^— -. (49) 

4 s ^(2»n!) 2 i u-u q -Y, u {q,u) K > 

Thus, one can write 

oo 

Z u (p,uj) = e~P 2 / 2 J2<j n (uj,u)p 2n , (50) 
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with 



a n (uj,u) 



q 2n+l e -q 2 /2 

n dq uj ~ uj q - E u (q,u)' 



(51) 



4(2"n!) 2 

Using such expansion allows one to promptly get a numerical solution by iterations (see Fig. [4j) . The convergence of 




0.2 0.4 0.6 0.8 1 
Figure 4: Real part of the renormalized spin stiffness as a function of the disorder strength u in units of es/4. 

the iterative solution is rather good up to some value of the disorder concentration u c . However, when u — > u c , we 
find that d u ai(0,u) — » oo. Therefore, it would be desirable to derive an analytical solution in the neighborhood of 
u c . For convenience, we omit the arguments of a n in our notation in the next part. In general, 



du u 4 



2 p "P 2 /2 



dq qc~ q2/2 7 o(gP) ~ 1 dY> u (q,uj) 
' ' (w - uj q - E u (q, to)) 2 du 



or equivalently 



dq 



q 2n+l e -q 2 /2 



da n _ <7n_ u 

du ~ u 4(2™n!) 2£ "7 ~* (u - u q - E u (q, lu)) 2 



dE u (q,u) 



(52) 



(53) 



Introducing for simplicity 



Eq. (|56| and has the form 



F n = 4 / ^ 2 (54) 

Jo [uj-uj q -T, u {q,u))\ 2 



one finds 



2^^«+fe— • ( 55 ) 



du u 4(2™n!) 2 ^ + du 
Introducing a matrix notation 

uF rn 

-\-n 



2 m + n + 2 m\n\ ' 



Eq. (|55j) reads 



fc=i 



, u cVfc _ 2 n n\a r , 
du u 



(56) 



(57) 



Its solution is found by computing the inverse matrix to 



da n 2 n \ - i k 



(58) 



k=l 



Substituting this result into 



d(det(i?) 2 ) = 2 / ddet(B) + ^ ddet ( B ) da n 

yields 



du ^— : ' da n du / ' 

n— 1 / 



d(det(S) 2 ) 



-k(v), 



(59) 



where 



k(u) =-2det(B) 



9det(g) 
du 



„ , /„x v 2 - adet(B) 2 fe -"fc! , , 
-2det(B) £ ^ ' — tTfcg-^, 



n.k— 1 



<9a n un! 
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with 



which leads to 



d -^l = det(S)Tr(V^ 
au \ au 



= -- det(B)Ti(B- 1 - I) 
u 



and 

ddet(B) 

where 



OR 



m+k 



,fc=l 



m 2«H-fc+2 m !/c! acr„ 



9F m 
<9er„ 



2e 



^2n+2m+l e -3g 2 /2 

i. d<Z [w - w g - £(g,u;)] 3 ' 



Suppose that det(B) — > when u — > u c . In this case 
B^ 1 dct(B) remains finite, as well as k(u). This suggests 
that 



dan 
du 



oc. 



(60) 



when u — > u c , since the other terms are finite. Moreover, 
if k(u) is a smooth function around u c , such that k(u c ) w 
k(mo) for some uo from the neighborhood of u c , then 
according to Eq. ([591 there holds 



dct[B(u)] = y/ k{u c ){u c — u) + 0(u c — u). 



(61) 



It follows then from Eq. ifBT]) that u c w 
dct[_B(uo)] 2 /Ac(Mo) as long as uq — > u c - However, the 
analysis of the infinite dimensional matrix i? and its de- 
terminant is quite complicated, which forces us to use 
an approximate solution, where we keep only the first 40 
terms in the expansion, thus reducing the dimension of 
the matrices to 40 x 40. In the absence of Zeeman split- 
ting (g — 0), for u> = 0, and uq = 0.238 one finds, setting 
( B = 1, that det[B(u )j = 0.0551776 and k(u q ) = 9.7945, 
which yields u c — 0.238311 in excellent agreement with 
the numerical solution. The approximation also allows to 
check the validity of Eq. flSfi]), which yields a[ = -9.384 
at the point uq — 0.238 (here the prime stands for the 
partial derivative with respect to u). On the other hand, 
the numerical solution for the two points uq — 0.238 and 
u\ = 0.23801 yields Acti/Am = —9.463, which agrees 
reasonably well with the previous result. The main dif- 
ference stems from the fact that uq = 0.238 is rather 
close to the critical point u c , where the derivative di- 
verges, so the value Au = 10 -5 is still rather large and, 
of course, computational errors and approximation with 
finite number of terms make the result not very precise. 
Furthermore, it follows from Eq. (|58|) that <7^det(i?) re- 
mains finite with u — > u c . Thus, 



cr n (u) - a„(u ) -- 
_ det[B{u )}a' n (uo) 



dv 



dcr n {v) 
dv 

dv(u c 



v) 



-1/2 



(62) 



, . , . , , , , det[B(u)l .„„. 
a n (u) = a n [u c ) + [a n {u ) - °"w( M c)J ^fg7^ y\ ( 63 ) 

after performing the integration, where <J n (u c )—a n (u ) = 
2{u c — Uo)(r' n (uo). From this analytic solution one may 
observe that a n and, consequently, S(p, ui) acquires an 
imaginary part when u > u c . In particular, considering 
n = 1, for the case at hand cti(wo) = —0.161742 and 
eri(w c ) = -0.167576. Defining 



a 



lim 



2a' 1 (u) det[B(u)] 



(64) 



the value of a can be evaluated without any fitting pa- 
rameters directly from Eqs. f58|) and l(60|) . which yields 
a = 0.331. It follows directly from the above that the 
renormalized spin stiffness now obeys 



p s (u) = p s {u c ) + ae By /u c - u, 



(65) 



where p s {u c ) = es(cri(u c ) + l/4). Both numerical and an- 
alytic results for p s {u) are plotted in Fig.0 which shows 
that the analytic solution remains in excellent agreement 
with the numerical one even for those values of u, which 
are far from the critical point u c . The whole behavior of 
the renormalized spin stiffness is very similar to the one 
obtained by Chalker et al.^- describing a spin glass phase 
transition. Moreover, such dependence of the renormal- 
ized spin stiffness as a square root function of a control 
parameter was already observed previously by Shender,— 
as well as by Avgin et al^ They considered the two- 
and tree- dimensional ±J Heisenberg spin glass model 
in a ferromagnetic ground state due to a strong external 
magnetic field. They found that for a certain value of 
the control parameter, p s (u) acquires an imaginary part. 
The real part of p s (u) is proportional to the spin wave 
stiffness, whereas the imaginary part is proportional to 
the damping of the spin wave excitations, thus signalling 
localization. It was argued that when the frequency of 
the spin-wave excitation u multiplied by its lifetime r 
is lot = Re[p s (u)]/lm[p s (u)] < 1, then the spin waves 
are completely localized. As we can see from the Fig.0 
the condition of localization is already satisfied for the 
values of the disorder strength starting from u = 0.3. 
The calculations presented in the Appendix D contain a 
strong indication that the Pauli susceptibility diverge at 
the point u = u c , suggesting a phase transition from a 
ferromagnetic ground state to a spin glass state,— since 
the spin waves become localized. 

Our discussion was mainly concerned with the static 
case lj — 0. However, our approach allows to find 
E(p, oj) for any given u. The dispersion spectrum in 
the self-consistent approximation then satisfies to — cj p — 
Re£(p,w) = 0. 
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Figure 5: (Color online) Real (black) and imaginary (blue) 
parts of the renormalized spin stiffness in units of e.b/4 as 
functions of the disorder strength u. (Square - analytic solu- 
tion, star - numerical) 



V. CONCLUSIONS & OUTLOOK 

This paper accounts for the presence of both disorder 
and interactions in a 2DEG at Landau level filling factor 
v = 1, whose ground state constitutes the well-known 
quantum Hall ferromagnet. The bosonization technique 
developed by Doretto et a/.— was employed in order to 
facilitate the treatment of both disorder and interactions 
in this strongly correlated system. The bosonization 
procedure consists in treating the spin wave (magnon) 
excitation as a boson such that the fermionic Hamilto- 
nian of the system can be approximately recast into a 
Hamiltonian expressed in terms of bosonic operators. As 
a consequence, the interaction between electrons up to 
RPA level was incorporated within the bare propaga- 
tor that represents the free boson. The intent was then 
to identify a disorder-driven quantum phase transition 
to a non-ferromagnetic state by analyzing the behavior 
of the renormalized spin stiffness as a function of the 
disorder strength, which itself corresponds to the ratio 
squared of the disorder energy scale to the Coulomb en- 
ergy one. To achieve this aim, firstly, we derived the 
bosonic expression for the Hamiltonian of the system. 
In the second stage, the focus was on seeking out the 
disorder self-averaged Green's function, which is the full 
bosonic Green's function averaged over the impurity po- 
sitions. Then, by using the Dyson's equation, we ob- 
tained a diagrammatic representation of the self-energy. 
The latter was subsequently computed within the frame- 
work of the low-density weak-scattering approximation. 
Low density means that the number of disorder atoms 
present in the system is taken to be much lower than the 
number of electrons, while the weak-scattering approx- 
imation signifies that the scattering potential induced 
by a given impurity atom is weak, such that only the 
first and second-order Born scatterings are accounted for. 
As a result, the self-energy corresponded to a single di- 
agram. Furthermore, the self-energy was evaluated in 
three different approximations: 1) the bare (long wave- 
length) approximation, which consists in using the bare 



bosonic propagator and keeping the lowest order terms 
in momenta, 2) the bare (full k) approximation, which 
uses as well the bare bosonic propagator but with all 
the momenta terms kept in the calculation and, finally, 
3) the self-consistent approximation, which uses the full 
disorder averaged Green's function instead of the bare 
one in the self-energy diagram. Then, the renormalized 
spin stiffness was determined by extracting the coeffi- 
cient of the quadratic term in the dispersion relation to- 
gether with the contribution from the self-energy. In the 
case of the bare (long wavelength) approximation, the 
spin stiffness was found to vanish linearly at the disor- 
der strength u p = 1. For the bare (full k) scheme, the 
spin stiffness also vanished linearly, but at the disorder 
strength u p = 0.7. These results suggest the occurrence 
of a disorder-driven quantum phase transition from the 
ferromagnetic phase to a paramagnetic one at the criti- 
cal value Up = 0.7. Lastly, the self-consistent calculation 
revealed a completely different behavior: the real part 
of the renormalized spin stiffness also initially decreases 
with increasing the disorder strength u, but then it satu- 
rates without reaching zero beyond a critical value u c , at 
which it (and the self-energy) acquires an imaginary com- 
ponent. According to the Shender criteriumjii the spin 
waves become completely localized when the imaginary 
part of the renormalized spin stiffness becomes larger 
than the real part, which occurs in our system for u > 0.3 
(see Fig.EJ). 

The physical mechanism behind a phase transition 
from the ferromagnetic ground state can be understood 
by considering electrons completely filling the lowest Lan- 
dau level [y = 1) in the presence of some inhomoge- 
neous electrostatic background (disorder). Then, for suf- 
ficiently strong impurity potential, by adjusting the elec- 
tron density to the electrostatic background, the system 
would gain more energy than is needed to rearrange the 
spin configuration. In this case the ferromagnetic state 
does not minimize the total energy of the system and a 
phase transition should take place. This quantum phase 
transition could be detected by calculating the behavior 
of the magnetic susceptibility as a function of the disor- 
der strength. A sharp peak is anticipated at the transi- 
tion point. In particular, if the energy cost for exciting a 
spin wave is less than the gain in the electrostatic energy, 
then the renormalized spin stiffness becomes negative and 
the system undergoes a phase transition to a paramag- 
netic state with zero local magnetization. On the other 
hand, as it was argued by Rapsch et aL,— in the case of a 
smoothly varying impurity potential, keeping nonzero lo- 
cal magnetization is still energetically favorable and the 
electrostatic energy is lowered by screening the impurity 
potential due to the formation of spin textures. At strong 
disorder such phase would correspond to a spin glass and 
the spin textures might be considered as the localized 
spin waves. Thus, the character of the phase transition 
might depend on the nature of the disorder. The cal- 
culations performed within our model indicate that the 
Pauli susceptibility diverges at the same critical point 
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of the disorder strength u c , where an imaginary part of 
the renormalized spin stiffness appears, thus suggesting 
a phase transition to a spin glass phase. 

Our approach can be extended for the case of bilayer 
systems in the presence of disorder. In fact, Fertig and 
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Appendix A: THIRD ORDER DIAGRAM 

Let us now evaluate the 3 rd order diagram contribution to the self-energy. Its diagrammatic representation is shown 
in Fig. 6. Algebraically, we have 

£( 3 > (p, u) = N lmp J2 U e (d ~ P, p)G(q, W )(7 e (q' - q, q)G(q', u;)U e (p - q', q') 



q q' 

xG(q»[/(p-q')/(p-q',q'). 



(Al) 



Here again, the impurity potential is short-range E/(q) = constant = U/A, By replacing all the functions defined 



(3) 

Z(p) 




U e (q-p,p) 



U e (p-q',q') 



Figure 6: Diagrammatic representation of the third order diagrammatic contribution to the self-energy, 
previously, we find 



E< 3 >(p,u,)=JV imp ^ 

q,q' 

2ie 



2 ^ e -K(q-P)iV4 si n f (q-p)Ap 



G(q» 



fq '- q)|2/4 dnf (q/ ~ q)Aq 



G(q,w) 

2ie -K(P-q')l 2 /4 sill f (P-q / )Aq- 



x sin 



U 
A 

q'Aq 



Term a 



N imp J2 ( 7 ) (2») 3 e-l^-P)l 2 / 4 e-^P-q')' 2 / 4 e-'^'-^' 2 / 4 sin f<l P 
|G(q,w)G(q , ,w) 



sin 



pAq' 



(A2) 



Term b 



Terr 



It is clear that Terms a and b are symmetric under the interchange q <->■ q' while Term c is antisymmetric. Thus, one 
has 



S (3) (p,c) = 



(A3) 



This result holds true in both the full bare and self-consistent approximations. 

As a matter of fact, due to the antisymmetric property of the wedge product within the sine term, it turns out that 
all odd order terms vanish. 
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Appendix B: DETAILED DERIVATION OF THE SELF-ENERGY 

We prove here the expressions (|39| and (|40|) . We begin with the expression of the self-energy given by Eq. J32 



E(p, W ) = An lmp U 2 r^g e -W 2 /2 e -W 2 /2 G(?)W) /'^e^^ cos ^[l-cos(^sin0)]. 



One first deals with the polar integral, 

fj g e^ c °^[l - cos^gpsin*)] = \ fj g 
The two terms are evaluated separately. For the first term, one must note that^ 

oo 

e £Vcos0 = h (f qp) + 2 £ I n (£ 2 qp) cos(n0), 

n=l 

such that 



e £ 2 gP co S _ 1 g e ^ gP cos0 cos(£2(?psin( 



(Bl) 



(B2) 



qp cos 



2 ./o 2vr 



— cos(n^) = -I (£ 2 qp). 



(B3) 



The second term 



1 



2n 



1 



^ e £ ^ cos *cos(£Vsm^, , , , „ . 

2 ,/q Z7T 2 Jq 27T 2 



2tt 



— exp[£ 2 (7p • exp(i<; 



(B4) 



Substituting Eqs. (|B3|) and (|B4|) back into Eq. (|B2[) then yields the simpler expression, 

I-2TT 11 



1, 



1 



_r e t qP cos<p rj _ cos( ^ psm0) ] = [J (£» - 1]. 
/q Z7T Z Z 

Now, substituting the above back into Eq. IjBip we find 

S(p, W ) = 4n imp U 2 J™ g 9 e-^ 2 / 2 e-^ 2 / 2 G(g,c)i[/ (^^) - 1]. (B5) 

Appendix C: BARE (LONG WAVELENGTH) APPROXIMATION 

To evaluate the self-energy within the long- wavelength approximation, we must return to Eq. f29|) . Firstly, one 
remarks that the sine squared term in Eq. f29|) greatly simplifies, 



sin2 (qA£) w (^P) 2 = 1 [/a , . (q x p)]2 = 1 ( , 2|q x p|)2 = l, 4|q|2|p 
Then, substituting Eq. dCTJ) into Eq. (29]) yields: 



2 sin 2 , 



E(p,w) = n imp U 2 i 4 / ^ g 3 e -^> 2 /y e -^) 2 / 2 G (q,^) 
o 2tt 



The polar integral then turns out to match^ 3 - 



2 ^d> 
o 2tt 



e < pgco S ^ sin 2 | 



(CI) 



(C2) 



The series expansion for Eq. <|C3|) gives (a; = £ pg) 
I (x) - 7 2 (x) 



2?r 



sin z cf> = -[I (e 2 pq) - I 2 (£'pq)}. (C3) 



u 


x 2 




\x 2 









(C4) 
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Eq. (|C2j) already holds a q 3 p 2 term and therefore a q 5 p 4 
term is not needed in the long wavelength approximation. 
Thus, one assumes that 



W 2 pq) ~ h(l 2 pq) w I- 



(C5) 



Moreover, the momenta are rescaled as q, p — > qj£,p/t. 
As a result, Eq. (|C2[i simplifies to 



l^{p,UJ) = ri P ' ' 



-p e 



dq q e 



3 p -9 2 /2. 



1 



UJ — UJ q + IT) 



where we replaced G°(q, uj) by its definition (see 
Eq. ([2T|) ). One can then make use of the identity^ 



1 



1 



x + irj x 



V in8(x), 



(C6) 



where V symbolizes the Cauchy principal value of the 
integral. Consequently, one has 

ReE(p, W ) = (H) up 2 e-^ 2 Vj o d q q -^— , (C7) 



— UJ„ 



ImS(p,w) = - f!£) 2 wp 2 e -P 2 /2 /"^ ^-9- .2 



X 7rS(L0 — UJq). 



(C8) 



Let us first examine the real part of the self-energy, 
which actually denotes the physical self-energy. 

It has been shown 1 — that in the long wavelength ap- 
proximation the bosonic dispersion relation for q can be 



written as 



e B 2 

^ q = .9 + -jq ■ 



(C9) 



The physical self-energy then becomes 
](^) = (^)\p 2 e- p2/2 7> [°°dq 



Re£( 



q 3 e -1 2 /2 



~u-g-^fq 2 ' 



4w _ Ag 

uj = — and g = — , 



(CIO) 

Let us then work temporarily with the new quantities 

(Cll) 

such that the self-energy is re-written as 

poo 

ReS(p,w) = e -fu p 2 e-P 2 / 2 V / dq- 
4 Jo L 



q 3 e- q > 2 



Now, one performs a change of variable in the q momen- 
tum: q — ► q = q 2 . One must note that qdq — d(q 2 )/2 and 
that the integration limits are not altered. Consequently, 
one gets 



ReS(p,w) = ^- U p 2 e- p2 / 2 V 



2-n f°° d( ? & 



-9/2 



o 2 uj- g-q 



(C12) 



A further change of the integration variable is performed 
q ^ k = uj — g — q, leading to 



ReSfoS) = pSe-P 3 /^ H^^-g- k) 

4 J Ld~g 1 



-(Q-g-k)/2± 

k 



£B 2 -p 2 /2 

— it pe p ' 



-(i-S)/2 



(C13) 



Term a 



-lv dkk^)e-^/ 2 



Term b 



Term a corresponds to the definition of the exponential 
integral function^ 



whereas Term b can be straightforwardly integrated, 



W-fl gfc/2 



■p j dk k— = I dfc e (*+a-S)/a 



Ei 



uj — g 



p r\(f ,a " 2 ' 



(C14) 



(C15) 
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where the shift of variable k k = k — (ui — g) was used 
in the first step. 

Thus, the physical self-energy becomes 



ReS(p,w) = ^upV p2 / 2 [-1 



(C16) 



The renormalized energy of the bosons is obtained by 
looking at the poles of the full disorder self-averaged 
Green's function, 



j p — ReS(p, uj) = 0. 



(C17) 



Consequently, in the long wavelength approximation, the 
renormalized dispersion relation takes the form 



to — g = p + u p e 
'ui-g 



V/2 r_! 



Ei [ <tL=J. \ e -(0-S)/3 



(C18) 



It is straightforward to notice that the renormalized 
spin stiffness, which corresponds to the coefficient of the 
p 2 term, is given by 



(C19) 



We now turn to the imaginary part of the self-energy 
given by Eq, l|C8|l , In the long wavelength approximation, 
the Dirac delta function becomes 



S(to (g + -^-Q 2 )) = ~9~ 9 2 )- 

(C20) 

By performing a change of variable in the q momentum, 
q — > q = q 2 and replacing Eq. (|C20|1 into Eq. l|C8j) one 
gets 

Im£(p,u>) = ~u p 2 e- p2 / 2 {uj - g)e- 2 ^- 9) ' tB . (C21) 

Finally, the scattering time, which amounts to the life- 
time of the bosonic excitation, is given by 



— =ttu p 2 e- p2 / 2 (u - g)e- 2 ^-^ €B . 



It is clear that t p — > oo when to — > g, i.e. low energy 
quasiparticles are long-lived, with finite lifetime induced 
by disorder. 



J 



Appendix D: PAULI SUSCEPTIBILITY 

The Pauli susceptibility in case of linear response is 
given by the Kubo formula: 

X«(x,x';t-t') =i(TA(x,t)S z (x',j')). (Dl) 
Using the Fourier transformation 



^ 2 (x,<)-^^(q,i)e lq - : 



(D2) 



the susceptibility can be written as 

X z M^-t-t') = i(T t SM,t)Sz(<{X))- (D3) 

On the other hand, the operators S z (q, t) can be written 
in the bosonized form^L 



SM, t) = - e-« 2 /4 J2 cos (fL^P) bl +p (t)b p (t), 

p 

where b p (t) — e lHt b p e~ lHt . Thus, after substitution 

xM;t-t') = ie - 2 / 2 E-(^)cos(^P: 

p,p' 

x (Tt6j 1+p (t)6 P (t)6^ +p ,(t , )V(^))- 
Evaluation of the expectation value yields 

Xzz(<l, q ; t - t ') = -ie q 1 2 2^ cos ( ] cos ' 



xG(p' + q', p; t - t')G(p + q, p'; t' - t), 



using the notation defined earlier in Eq. 1(13)1 . Expand- 
ing the Green's function G(p,q;i — t') into the Born 
series and performing the disorder averaging one recov- 
ers the translational invariance (Xzz(q.,<l';t ~ £'))im P = 
^q+q',oXzz(q, t — t'). Moreover, performing the Fourier 
transformation in the time variable t and introducing 



(C22) P(p,q;u,,e) = ^cos 



qAp' 



x (G(p' - q,p;w + e)G(p + q,p';w)) imp 
the susceptibility is 



Xzz(q, e) 



iAe-^l 2 
(2tt) 3 



dp/ dwcos (^y^) P(p,q;w,e). 



(D4) 



17 



In the self consistent approximation the function P(p, p'; u>, e) obeys^ 

'pAp'\ AN; 



P(P, p'; w, e) = G(p, w + e)G(p + p', w) 



cos 



-^T J dqJ7 e (p-q,q)C/ e (q-p,p + p') J P(q,p';w,e) 



2 y ' (27r) 5 

We are interested mostly in the static susceptibility \ = lim c _o Xzz(0, e)- Thus, in particular 



P(p,0;^,0) = G 2 (p,^) 
A spherically symmetric solution satisfies 

Let us introduce a new function 



1 



¥ imp 



(2tt)s 



J dq[/ e (p-q,q)£/ e (q-p,p)P(q,0;w,0) 



,2 ,-p 2 /2 



dg qe- q2 / 2 (I (qp)-l)P(q,0;u,0) 



H(p,u>) = P(j>,0;u,0)G- 2 (p,ujy, 



then Eq. i|D5|) can be rewritten as 



ff(p,u;) = 1 + -e%e-P 2 / 2 J^ dq qe-" 2 / 2 (I (qp) - l)G 2 (q,ui)H(q,ui), 



or explicitly 



/•OO 

H(p,o;) = l + |e|e-f 2 / 2 y dq 



-„ 2 /2 ^o(gp)-i 

[w - w g - S„(g,o;)] 



2 /%,u;). 



(D5) 



(D6) 



(D7) 



(D8) 



Notice that Eq. 1|D8|) has the same form as Eq. lj52]) but with w) instead of d u T, u (p 7 w), which is known to diverge 
d u T, u (p, 0) — * oo when u — ► u c . In the next part we will demonstrate that H(p, 0) also diverges, 0) — > oo when 
14 — > u c . 



We are looking for a solution in the form 

oo 



r 



and 



71=1 



(D9) 



Substitution of Eq. (|D9|1 into Eq. I|D8|) yields an expres- 
sion, which looks similar to the equation previously ob- 
tained (see Eq. (|55|) ). 



Equivalently 



fe=0 



(-1)* 
2 k kl 



K, 



fc+n ■ 



h n (uo) 



uK n 



4(2™n!) 2 4(2«n!) 



fe=i 



2 n + 2 n\ ' 



(D13) 



(D14) 



fe=i 



where , ra was defined in Eq. (|56|) . The solution is 
where the function F n was defined earlier by Eq. (ED and found by computing the inverse matrix to Eq . QgQ and 



Notice that 



dq 



q 2n+l e -q 2 /2 



has the form 



[w-w q - E u (q,u>) 



fe=0 



(Dll) 



(D12) 



2" n it 



4n! ^~™> fc 2 fe fc!' 

fe=i 



(D15) 



Therefore, 



2~ n u ^ „_ x K k 



P{p, 0; w , 0) = G 2 {p, U )[l + e-^' 2 E^^E 



fe=i 



: 2 fc /c! 



(D16) 
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and 

/>oo roc 00 roc c\ — n 00 / ' 

jf P(p,0;^0)pdp = j o G\p,u;)pdp + J2l V n+1 G?(p, W )dp^ g S" 1 ^, (D17) 
if the integral is convergent. Otherwise, it has to be regularized, which we won't consider here. This leads to 

n— 1 n.k— 1 

which can be further simplified by means of some algebraic transformations, 

l°°P( P ^,0)p d p = F 0+ ± (^< fe -« + ^-^ + |^-^). (D19) 

n,fc=l 

I 



Despite the simplifications, the above expression is diffi- 
cult to evaluate analytically, as well as numerically. How- 
ever, since most of the terms there involve the inverse ma- 
trix, it is reasonable to suppose that if u> = it diverges 
with u — > u c as 

/>oo 

/ P(p,0;0,0)pdp ~ det[-B(u)] _1 . (D20) 
Jo 

On the other hand the susceptibility is given by 

iA f°° f°° 

x(e) = _ (2^)2 y dw j p (p>°^>t)p d p- (° 21 ) 



Thus, considering e = 0, we see that the integrand is 
divergent at u> = with u —* u c , which is definitely not 
enough to infer the divergence of the integral itself, but 
can be considered as an indication to such possibility. 
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